Inelastic Quantum Transport and Peierls-like Mechanism in Carbon Nanotubes 
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We report on a theoretical study of inelastic quantum transport in (3m, 0) carbon nanotubes. By 
using a many-body description of the electron-phonon interaction in Fock space, a novel mechanism 
involving optical phonon emission (absorption) is shown to induce an unprecedented energy gap 
opening at half the phonon energy, huJo/2, above (below) the charge neutrality point. This mecha- 
nism, which is prevented by Pauli blocking at low bias voltages, is activated at bias voltages in the 
order of hu)o . 
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Since their discovery, carbon nanotubes (CNTs) have 

attracted much attention due to their outstanding me- 

i—i - 

chanical and electrical properties [lj . Depending on their 
helicity, CNTs exhibit either semiconducting or metallic 
behavior and for low resistance contacts, ballistic trans- 
port is observed in the low bias regime. At higher bias, 
several experimental works report on current saturation, 
attributed to inelastic (inter-band) backscattering by op- 
tical phonons Q{. Electron-phonon (e-ph) interaction 
thus severely limits the exceptional properties of CNTs 
as ballistic conductors. Aimed to the explanation of the 
experimental data 0- 0- 0] ■ several theoretical studies of 
electronic transport in CNTs inpresence of e-ph inter- 
action were performed H, IE S S S [XZ3] • They include 
simulations based on the Boltzmann equation 0, E3 , the 
Fermi Golden Rule |(j and the use of a diagonal self- 
consistent Born approximation Q. Besides, in order to 
explain the quantitative disagreement between theoreti- 
cal and experimental estimations of inelastic mean free 
paths, a hot phonons scenario was recently proposed 0. 

On the other hand, electron-phonon coupling in low- 
dimensional systems can lead to more subtle effects and 
important corrections to both phonon and electronic 
bandstructures. Fundamental examples are the Kohn 
anomaly and the Peierls transition [l2L fl3| . Con- 
nected to the latter, and by using density functional 
theory, Dubay and coworkers found a softening of the 
mode with A\ symmetry in metallic CNTs, not suf- 
ficient however, to produce a static lattice distortion 
|l4j . Notwithstanding, as explored within the linear re- 
sponse regime, the corresponding phonon-induced (time- 
dependent) electronic bandstructure changes result in 
strong modifications of the Kubo conductance 0. This 
points towards the importance for a quantum mechani- 
cal treatment of inelastic transport whenever high energy 
(optical) modes are activated. 

In this Letter, by using a full quantum description of 
the joined processes of tunneling and phonon-assisted 
transport, the coupling between electrons and optical 
Ai(L) phonons in (N = 3m, 0) zig-zag tubes is shown 
to result in an energy-gap opening at fiioo/2 above (be- 
low) the charge neutrality point (CNP), owing to phonon 
emission (absorption). This novel many-body mechanism 



is activated when driving the system out of equilibrium. 

For simplicity, we consider an infinite CNT and al- 
low the electrons to interact with phonons only in a 
central part of length L. The Hamiltonian of the sys- 
tem is written as a sum of an electronic part, a phonon 
part and the e-ph interaction term. The electronic 
part is described by an effective 7r-orbitals model H e = 
~'yoJ2{i. 3 ) [ c t c j + h.c], where c\ and c t are the cre- 
ation and annihilation operators for electrons at site i, 
7o = 2.77eV is the transfer integral restricted to near- 
est neighbors 7r-orbitals 1]. The phonon term is given 
by H p h — hujQb + b, with b + and b the phonon oper- 
ators. The remaining contribution is given by a Su- 
Schreiffer-Heeger Hamiltonian [l^, describing the e-ph 
coupling to a non-local vibrational eigenmode. This is 
found by assuming a phonon modulation of the hopping 
matrix elements |l5l 1 1 6| , keeping only the linear correc- 
tions to the atomic displacements from equilibrium, i.e. 
7i,j = 7o + a ^i,j ' &Qi,ji where <5jj is a unit vector in the 
bond direction, 8Qi j is the relative displacement of the 
neighboring carbon atoms, and a is the e-ph coupling 
strength defined as the derivative of 70 with respect to 
the change in the bond length. Further quantization of 
the atomic displacements gives the e-ph interaction term: 
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where the e-ph matrix elements can be written as 7, 

a.y/h/ (2mujo)Si - j ■ {ei — ej), e, is the phonon mode eigen- 
vector which gives the atomic displacements. The sum- 
mation in the right hand side corresponds to nearest 
neighbors within the vibrating region of the CNT. For 
the phonon mode considered here (Ai(L)), A and B-type 
atoms move out of phase in the direction parallel to the 
tube axis (see scheme in Fig. 

To compute inelastic quantum transport, we use the 
approach introduced in Refs. 0,0], which has been ap- 
plied to a variety of problems including vibration-assisted 
tunneling in STM experiments [19(, transport through 
molecules [2iJ, and resonant tunneling in double bar- 
rier heterostructures 1231 . This scheme allows an exact 
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non-perturbative treatment of the problem of one elec- 
tron transport in the presence of e-ph interaction. This 
is achieved by an exact mapping of the many-body prob- 
lem into a single particle problem in a higher dimensional 
space (the e-ph Fock space), each considered phonon 
mode adding one extra dimension to the problem. The 
key idea behind this approach can be visualized after re- 
writing the interacting Hamiltonian in an appropriate ba- 
sis for the e-ph Fock space (one electron plus phonons). 
In this equivalent multichannel one-body problem, the 
asymptotic states (in the non-interacting leads) include 
both the electronic and the vibrational degrees of freedom 
and the total energy is fully conserved. When consider- 
ing a single phonon mode, these asymptotic channels can 
be labeled by using two indexes: (X, n), being X = L,R 
the index associated to the corresponding electrode (left 
or right) and n the number of phonon excitations in the 
system. Accordingly, the transmission ^pr n )— >tKm) and 
reflection probabilities R^x n)-*(Y m) between the differ- 
ent channels are computed by using standard Green's 
functions techniques [22]. Using these probabilities as 
inputs, the non-equilibrium electron distributions in the 
leads (at finite temperature and bias voltage) are eval- 
uated self-consistently using the procedure developed in 
Ref. |2l|. The electronic current through the device can 
be obtained from these self-consistent distributions which 
take into account the Pauli exclusion principle for the dif- 
ferent competing processes. 

Let us first explicitly establish how the Fock states 
are connected by the Hamiltonian and simplify the prob- 
lem by using a mode decomposition. The matrix ele- 
ments of the e-ph Hamiltonian are essentially given by 
the projection of the bond direction on the relative dis- 
placements of the atoms from their equilibrium positions, 
' (ci — Ej). It can be easily proven that the e-ph ma- 
trix elements 7?"? have values 7q~ p/i = —2a>y/h/ (2mcJo) 
for bonds that are parallel to the tube axis, and 

— lo' ph cos(7r/3) for bonds that are tilted with re- 

spect to the tube axis. 

For the phonon mode considered here, instead of solv- 
ing the Hamiltonian in real space, one uses a mode space 
approach [3, |23j ■ The idea is to resort to a unitary trans- 
formation that diagonalizes the electronic Hamiltonian 
for each layer perpendicular to the tube axis. The re- 
sulting eigenstates will correspond to different circum- 
ferential modes that can be used to build an alterna- 
tive basis for the description of transport through the 
nanotube. In absence of static disorder, the electronic 
Hamiltonian H e does not couple the different modes, 
and the different subbands correspond to linear chains 
with alternating hoppings 70 and 7 9 = 2jq cos(qn / N) 
(q = 0, l...,N — 1) with dispersion relations: £^°'(k) = 
iy^Q + 7^ + 2707^ cos(3fca cc /2). When N is an inte- 
ger multiple of three, the metallic subbands correspond 
to q = N/3, 2N/3. Furthermore, since the e-ph Hamil- 
tonian H e _ p h does not couple different circumferential 
modes neither, one is left with N independent problems. 
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FIG. 1: Total transmission probability as a function of the 
energy of the incident electrons for no — 0. Inset: Same in- 
formation for no = 1. The results correspond to a (24, 0) tube 
where interaction with an Ai(L) mode inducing displacements 
along the axis direction (drawing) is included. Note that a 
total energy e = huio/2 corresponds to an electronic kinetic 
energy E — hwo/2 for n = and —huio/2 for n = 1. 



Each of them corresponds to a one dimensional mode 
space lattice that, when coupled to phonons, gives a two- 
dimensional problem in Fock space. 

In Fig. ^ one reports the total transmission proba- 
bilityT(S) = J2 n T (L (E) as a function of the 

incident's electron energy E, when there are initially no 
phonon excitations in the system, uq — 0. All the ener- 
gies are expressed in units of 70, the value of Hujq is taken 
as O.O770 |1J] and a = oto ~ 7eV/A is estimated from 
|24j . For short nanotube lengths [L = lOnm), the main 
feature is the occurence of a transmission dip centered 
at E = Iuuq/2 above the CNP, that progressively deep- 
ens with CNT length, to reach a full gap at L = lOOnm, 
and with a width of approximately 3jQ~ ph . Apart from 
this fundamental gap, the effect of e-ph interaction re- 
mains small. The transmission remains mostly elastic in 
all the energy range shown in the plot. However, the re- 
duction of the transmission probability in the gap region 
corresponds to a complementary increase in the inelastic 
backs cattering by phonon emission. 

Figure ^inset) shows the transmission probability in 
the situation where one phonon is already available for 
scattering before any additional charge enters into the 
sample. In this case, in addition to the expected en- 
largement of the gap at E = hujo/2 by a factor y2 due 
to stimulated phonon emission, another gap develops at 
E = —Hujq/2. This transmission gap is complemented 
by an increase in the inelastic backs cattering by phonon 
absorption. 

To provide a clear physical picture of such phenomena, 
let us consider that the sample, i.e. the region where the 
e-ph interaction is allowed, extends to infinity. The Fock 
space for the coupled e-ph system can then be expanded 
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FIG. 2: Scheme for the unperturbed (H e ~ p h = 0) dispersion 
relations corresponding to n = 0, 1, 2 for the mode q — N/3. 
The states marked with identical symbols correspond to de- 
generate states of the coupled e-ph system. These degenera- 
cies are lifted by the e-ph interaction leading to the opening 
of energy gaps at hwo/2 above (below) the band center (see 
right inset). Upper inset: Representation of the many-body 
Hamiltonian in Fock space for a circumferential mode q (note 
the two-layer periodicity). Solid circles represent states in 
Fock space while the lines are off-diagonal matrix elements. 



in terms of the basis states {\l q ,n) = \l q ) ® \n)}, where 
\l q ) is the q-th circumferential mode localized in the I— th 
layer and \n) corresponds to the state with n phonons in 
the system. Close to the CNP, only two circumferential 
modes participate in conduction (q = N/3,2N/3). The 
e-ph Hamiltonian for one of those modes is represented 
in Fig. |21 (upper inset). Alternatively, the Fock space 
for an individual mode can be expanded by using plane 
waves instead of localized states for the electronic part. 
This leads to basis states of the form \k,n) = \k) <g) \n) 
where |fc) is a plane wave in mode space with wave- vector 
k along the axis direction. 

Let us consider the circumferential mode with q = 
N/3, though a similar reasoning holds for q = 2 A/3. In 
absence of e-ph coupling, one gets disconnected chains 
associated to the different values of n, each one with 
hoppings 7 g = 70 and unperturbed dispersion relations 

£1? = £n \k) (see Fig. |2J). Of particular interest are 
the states \k + — n/2a + \Sk\ ,0) and |fc + — K, 1) (black 
circles), where a = 3a cc /4 and K = n/a. They have 
the same total energy e^\k + ) = e^(k + — K) = fko a /2. 
When the e-ph interaction is switched on, the lattice pe- 
riod of the many-body Hamiltonian is doubled because 
of the spatial periodicy of H e _ p h. The new lattice vector 
is K = tx I a instead of 2ir/a. The mentioned Fock states 
are now mixed by H e _ p h, i.e. (k + — K, 1| H e -ph |^+>0) — 
AY~ ph = (3/2) 7 o" p/t /0@. The degeneracy is thus 
lifted, giving rise to the opening of an energy gap in £o(k) 
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FIG. 3: Total transmission probability at E — huo/2 as a 
function of the tube length for no — 0, the other parameters 
are the same as in Fig. The inset shows the saturation 
value of the total transmission as a function of a 2 . 



and £\{k) of width 2 |A7 e " p/l | around e = Hujo/2 (right 
inset in Fig. [2J. Accordingly, an incoming electron with 
a wave- vector k in the Fock state \k,n) will contain, as 
k approaches k + (or k + — K), an increasing admixture 
of \k — K, n + 1) (or \k + K, n — 1)) leading to a Bragg- 
type of inelastic scattering. This mechanism is due to a 
modification of the translational symmetry of the system 
driven by the e-ph interaction, with no static distortion 
of the lattice. A similar effect holds for the other states 
marked as symbols of the same kind in Fig. and for 
the sub-band with q = 2 N/3. This results in the trans- 
mission gaps shown in Fig. ^ which are in quantitative 
agreement with the prior analysis. 

Figure[3]shows the transmission minimum at the center 
of the phonon emission gap T(E — Hujq/2) as a function 
of the nanotube length. Different curves correspond to 
different values of the parameter a entering into the e- 
ph matrix elements li~J h ■ From a simple argument, one 
can predict the existence of a region where the transmis- 
sion decays exponentially with the tube length (tunnel- 
ing through the gap), i.e. T oc exp(— £/£), with a decay 
length £ that is inversely proportional to the energy gap 
A-f e ~ ph . Another interesting feature is the observed sat- 
uration for longer tubes. The saturation value T sat as a 
function of a 2 (inset of Fig. is associated to the small 
inelastic component of the transmission. 

It should be noted that the opening of the transmission 
gap will be prevented by Pauli blocking at low bias volt- 
age. In order to activate it, the system has to be driven 
out of equilibrium by applying a sufficiently high bias be- 
tween the voltage probes. To explore the consequences 
of this phenomenon in the current-voltage characteristics, 
the effect of the bias voltage is introduced on H e by mod- 
ifying the onsite energies of the 7r-orbitals. The poten- 
tial drop is assumed to be equally distributed at the two 
contacts. The current-voltage curves calculated from this 
simple model for a = ao and a = 2ao and zero tempera- 
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V [V] 

FIG. 4: Current vs. bias voltage for a (24, 0) 150 nm long 
tube for a — cto (main frame) and a = 2ao (upper inset). 
Solid lines correspond to a fixed phonon population no = 
0. The dashed, dash-dotted and dotted lines correspond to 
the effective thermal distributions described in the text. The 
mean number of phonons corresponding to these curves are 
shown in the lower inset. 

ture are shown in Fig. ^ (solid lines in the main frame and 
inset respectively). The main feature that results from 
the gap opening is the onset of a current plateau observed 
at V ~ frujQ whose width scales linearly with a. The sub- 
sequent linear increase in the current is due to the lack 
of other ingredients in our model such as electronic cou- 
pling with other phonon modes (such as those producing 
inter-band backscattering) or electrostatics effects. Moti- 



vated by recent studies that suggest the possible build up 
of a strong out-of equilibrium phonon population at high 
bias voltages 0, IE El > a natural question is how such 
scenario would modify our results. The answer might 
be provided by a self-consistent scheme to determine the 
out-of equilibrium phonon population in the system, but 
this is beyond the scope of the present study. Instead, 
a thermalized phonon population with an effective tem- 
perature ksT e ff is assumed and kpT e jf{V) is supposed 
to increase linearly for eV > Hloq |3. ItiTlCj. The depen- 
dence of the mean numbers of phonons with the bias volt- 
age for different slopes in k B T eff {V) (lOkK/V, 20kK/V 
and 30kK/V) are shown in the lower inset. The corre- 
sponding currents for a 150nm long tube are shown with 
dashed, dash-dotted and dotted lines for a — ao (Fig. 
QJ and a = 2ao (upper inset). A substantial decrease 
in the currents when the phonon population increases is 
observed. Furthermore, the dotted lines reveal the ap- 
pearance of a region of negative differential resistance. 

Semiconducting zig-zag CNTs with a bandgap smaller 
than hujQ are also found to display the phonon-induced 
transmission gap (not shown here). On the other hand, 
for CNTs other than zig-zag, this phenomenon will be 
driven by the relevant optical phonon mode producing an 
instantaneous dimerization pattern of appropriate period 
in the axis direction. In contrast to the Peierls transi- 
tion, where a mean field description is feasible ^5|j the 
proposed mechanism does not involve a static lattice dis- 
tortion and is activated at high bias voltage. 

This work was supported by the French Ministry of Re- 
search under grant RTB: PostCMOS moleculaire 200mm. 
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